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We study order parameter fluctuation effects in the superconducting state as a possible precursor 
to the pseudogap phenomena. Using a low-energy effective theory in the d-wave BCS model, we 
self-consistently calculate the single-particle properties. We find that the fluctuations reduce the 
spectral gap, and the shape of the gap is deformed, showing reduced slope of the gap at the nodes. 
We also show the angle-dependence of the quasiparticle lifetime due to the fluctuations. 



Introduction. — Despite much experimental evidence of 
pseudogap phenomena in the underdoped cuprates, their 
microscopic mechanism is not understood.!!! However, a 
pairing precursor as the origin of the pseudogap is one 
prominent possibility.B An active endeavor has been to 
incorporate strong pairing fluctuations to account for the 
pseudogap phenomena, especially above the supercon- 
ducting critical temperature T c .cla In this report, we take 
a different tack and study the order parameter fluctua- 
tion effects in the superconducting state. From a phe- 
nomenological standpoint, the pseudogap state can be 
considered as a superconductor whose phase coherence is 
destroyecLby strong phase fluctuations whereas the gap 
is robust .13 Therefore, it will bare similarity to a super- 
conducting state with strong order parameter phase and 
amplitude fluctuations. Here we do not attempt to repro- 
duce pseudogap phenomena since we study fluctuations 
in the weak-coupling BCS theory below T c and do not 
consider vortex pair unbinding transition, but much of 
the qualitative trend is expected to pertain to the pseu- 
dogap state. Using the effective low-energy theory ap- 
proach, we may describe the problem with relatively few 
physical parameters and separate the effect of order pa- 
rameter phase and amplitude fluctuation effects oQ 

We address the following two issues concerning the 
fluctuation effects. First, we examine the effect of the or- 
der parameter fluctuations on the size of the spectral gap 
below T c . It has been shown that the fluctuations reduce 
the magnitudje-|£)f the order parameter and the critical 
temperatureJSBo In this paper we also show the reduc- 
tion in the spectral gap in a wide temperature range. For 
simplicity we do not the include Coulomb interaction or 
disorder although they would alter the form of the order 
parameter fluctuations and further modify our result if 
properly included.ua Secondly, we examine the angular 
variation of the fluctuation effect. Angle-resolved pho- 
toemission spectroscopy (ARPES) data on underdoped 
cuprates show that above or below T c the shape of the 
gap near the TOde significantly deviates from the simple 
d-w&ve shapeJiS'EII We show that this may be due to the 
fluctuation of the phase rather than the amplitude of the 
order parameter, and moreover that the amplitude fluc- 
tuation effect is the strongest near the antinode. Also we 
discuss the angular variation of the quasiparticle lifetime. 
Low-energy effective theory. — We consider the weak- 



coupling mean-field BCS theory in which the pairing po- 
tential gives a d-wave order parameter which is effective 
only within the momentum thickness 2 A <C hp around 
the Fermi surface. Then the effective fcrmion Hilbert 
space is a thin momentum shell of a characteristic thick- 
ness 2A around the Fermi surface where A <§C kp. For 
convenience, we coarse-grain the momentum shell into 
small boxes and label them with an angular variable <fi. 
The effective action of the fermions and the order param- 
eter is 
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where \P is the superconducting order parameter intro- 
duced via Hubbard- Stratonovich transformation to de- 
couple the pairing interaction. In the above we assume 
a pairing potential of the form V{4>,4>') = g w{<f))w((j)') 
where w(<f>) — cos 2<f> and g < which produces a d- 
wave order parameter. It is understood that in writ- 
ing c CT (0;k, r), the momentum k lives only inside the 
small box labeled by the angular variable <j) near the 
Fermi surface. In order to explicitly separate the or- 
der parameter phase and amplitude degrees of freedom, 
we re-express \&(x, r) = A(x, T)e l6, ( x,r ) where A(x, r) 
takes a real value. In the mean-field approximation, 
we replace A(x, t) with Ao and obtain a d-wave gap 
A(^) = Ao cos 2<f> using the following self-consistent gap 
equation: 
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where £k = k 2 /2m — fx. The momentum summation 
above is constrained by the condition |£k| < vfA. Here 
we consider the fluctuation around the mean-field value 
and re-express A(x, t) = A + d(x, r). Then we perform 
a gauge transformation ^(x, r) = c CT (x, r)e _,6, ^ X:T ^ 2 , to 
couple the phase fields to the fermions explicitly. 
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The resulting effective action is expressed in terms of 



the Nambu spinor notation, ip 
Si, with 
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which amounts to summing over the ring diagrams in cal- 
culating the order parameter correlation function. The 
fermion self-energy can be obtained self-consistently from 
the Dyson equation. 

We first evaluate the correlation functions of the am- 
plitude fluctuations: 
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and d(is, q) = a x d(y, q), with Go the bare Green's func- 
tion for the neutral fermions. In the above, we approxi- 
mate £k = k 2 /2m — fi as £k ~ w F(|k — kp) and vp{(f>) is 
the Fermi velocity in the <j) direction. 

In building this effective theory we have not consid- 
ered vortex pair unbinding which leads to the Kosterlitz- 
Thouless transition. This is justified well outside the fluc- 
tuation regime. We also assume that we are in the tem- 
perature range where the BCS mean-field theory is justi- 
fied, namely, that ST /JL Ao /Ep in a two dimensional 
clean superconductor.E 2 ] One thing we observe from the 
form of the effective theory is that the strength of the 
coupling between fermions and the amplitude fluctua- 
tions has an angle-dependence; the amplitude fluctuation 
effect is suppressed near the gap nodes, as is evident in 
Eq. (|ll| ) where the self-energy correction is multiplied 
by a factor of w 2 (4>) . The strength of the coupling to the 
phase fluctuations is not suppressed at the node, how- 
ever. 

In studying the finite temperature superconducting to 
normal state transition, it should be sufficient to consider 
only the static fluctuations of the phase, and we may 
suppress the time-dependence in 9 and d and retain only 
the spatial fluctuations. Eq. ([|) is then modified as 
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The simplified form above, however, does not produce 
reliable results near zero temperature. 

Fermion single-particle properties. — In evaluating the 
quasiparticle self-energy by perturbative expansions, we 
take advantage of the fact that the effective theory re- 
sides in the thin shell around the large Fermi surface and 
select the diagrams which are of leading order in A/kp, 



From Eq. (g) we obtain (d(q) d(-q)) rin g = T/(a + b q 2 ) 
where the temperature-dependent coefficients a and b can 
be evaluated by carefully expanding Tldd in q from Eq. 
(|J). If we only consider the spatial fluctuations in the 
order parameter, the (d 0) terms are zero in the Gaus- 
sian approximation. Also the phase fluctuation has the 
following well-known correlation function!] 
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where n s {T) is the superfluid density at temperature T. 

Now we can determine the quasiparticle self-energy 
correction using the self-consistent Dyson equation, ne- 
glecting the vertex corrections: 
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where G is the full fermion Green's function, given self- 
consistently by G _1 = Gq 1 — X. In general the self-energy 
has both a momentum and frequency dependence, but we 
focus on the behavior of the self-energy near the Fermi 
surface, assuming that it varies smoothly near the Fermi 
surface. Therefore we neglect the ^-dependence so that 
the only momentum dependence is through the angle <f> 
on the Fermi surface. Then we can approximately obtain 
the self-energy: 
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where Cj and A can be calculated self-consistently by 
G- 1 (&k,w) = Gq 1 fo;k,w) 



2 



From the self-energy obtained in Eq. (|Tl|), by an- 
alytically continuing the frequency iui — > lu + iij, we 
can calculate various single-particle properties such as 
density of states (DOS), spectral functions, and single- 
particle scattering rates. In this paper we focus on DOS 
since it is gauge-iwariant and measurable via the tunnel- 
ing sp£ctroscopyc3 or the momentum-integrated ARPES 
datallil. We are especially interested in the angle-resolved 
DOS: 
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as it gives information about the angular variation of the 
fluctuation effect. 

Throughout this paper, we set the relevant energy 
scales A « 10A (T = 0)/vp and E F « 5Aw F so that 
we are well in the BCS weak-coupling regime; these rel- 
ative energy scales give much stronger pairing strength 
than ordinary superconductors but significantly weaker 
than the cuprates. With the energy scales so chosen, 
we may estimate the regime of the validity of the mean- 
field theory. If we apply the Ginzburg criterion, namely, 
Ao(T)J 2 3> (d(x) rf(x)), which may be estimated from 
Eq. (R), the mean-field theory breaks down only near 
T/T c ~ 0.98. Therefore, the BCS framework is reliable 
in most of the temperature range that we consider. 

In Fig. [j] we show the total DOS. As the temperature 
increases, the DOS peak is widely smeared. Very close to 
T c , the DOS peak has almost disappeared and the spec- 
tral gap is only manifested by the depletion in the DOS 
around w = as compared to the normal state DOS. 
Figure || shows the DOS peak position as a function of 
temperature; in most of the temperature range, we can 
interpret the peak position roughly as the spectral gap. 
This figure shows that the spectral gap is reduced in a 
wide temperature range due to the order parameter fluc- 
tuations. Near T/T c rj 0.98, the DOS peak structure has 
almost disappeared, and the DOS maxima do not have a 
meaning as the spectral gap. Therefore, we need to esti- 
mate the size of the spectral gap from the width of the 
DOS depletion in this case. It is difficult to study the evo- 
lution of the spectral gap through T c in this framework 
because the mean-field theory breaks down sufficiently 
close to T c as discussed above. Also the result near zero 
temperature is not reliable due to the negligence of the 
time-dependence of the fluctuations. 

The angle-dependent DOS peak near T c is shown in 
Fig. |. At T < T c , the DOS peak contour follows the 
d-wave shape. As the temperature approaches T c , we 
observe that the DOS is widely smeared to low-energy 
states especially near the antinode (4> = 0). Figure || 
shows that the shape of the angle-resolved gap (DOS 
peak curve) is deformed from the original d-wave shape 
near the node. We argue that the downward bending of 
the DOS peak curve near the node (0 = 7r/4) is due to 
the phase fluctuations since the amplitude fluctuations 
alone do not cause the downward bending as illustrated 



in the same figure. We find that the angle-dependence of 
the spectral gap near the node is 
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and the slope of the gap near the node is reduced. It has 
similarity to the shape of the gap obtained by ARPES 
on underdoped BigS^CaC^Og+i (Bi2212) in the super- 
conducting statejld although its microscopic origin is not 
well understood. 

Figure || shows angular variation of the scattering rate 
due to the order parameter fluctuation. We observe that 
the maximum scattering rate occurs near the antinode 
and also that the rate decreases as one approaches the 
node. This variation is due to the angular dependence 
of the order parameter magnitude fluctuations. This fea- 
ture may contribute to the anisotropy of thepquasiparticle 
scattering rate in cuprate superconductors .E-3 

Discussions and Conclusions. — Here we discuss the 
qualitative effects of amplitude and phase fluctuations. 
In Fig. 3, we find that if we omit the phase fluctuation 
effect, the apparent gap is enhanced. This is because 
the amplitude fluctuations tend to increase the gap mag- 
nitude. This can be understood from the Fermi liquid 
reference frame as following: The fermion self-energy cor- 
rection due to the amplitude fluctuations can be roughly 
estimated as 

E(p, u) « - J dvd 2 qG{p - q, u - v) (A(q, i/)A(-q, -v)) 
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and therefore the effective spectral gap is enhanced by 
the fluctuations as |A cff | 2 = (A(x)A(x)) = |A | 2 + 
(SA{x)SA(x)). 

The effect of phase fluctuation can be considered as a 
Doppler shift in the fermionic spectrum by ■ v s where 
v s ~ V0/m. Due to thermally fluctuating superfluid 
velocity, the DOS near the gap is now shifted since the 
energy levels at the gap nodes are enhanced. As a result, 
more states would be occupied near the nodes, and hence 
decrease in the slope of the gap nodes as shown in Fig. 3. 
On including the vortex pair unbinding transition of the 
BKT type, which gives stronger phase fluctuation effects, 
we can obtain a Fermi arc-like phenomenon .Qila 

In order to obtain the correct effect on the size of the 
spectral gap, both the phase and amplitude fluctuations 
have to be self-consistently taken into account. The total 
effect is reduction in the spectral gap as shown in Fig. 2 
and 3. However, the spectral gap may not be equal to the 
order par.ameter magnitude especially if the fluctuation 
is strong ,Ej and hence more careful study is needed to 
separate these two quantities. 

The results presented in this report are equally well 
applicable to any unconventional superconducting sym- 
metry. In principle, any superconductor would have a 
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window of temperatures near T c where such fluctuations 
are visible, depending on the pairing strength and the 
superfluid density. In case of underdoped cuprates, how- 
ever, it would be essential to include the effect of vortex 
pair unbinding in the pseudogap state, due to the small 
superfluid density. The result of this report may never- 
theless pertain to its superconducting state. Indeed, a r£=. 
cent observation of the deformed gap shape in Bi2212,EU 
which is only observed in underdoped regime, may be 
related to the order parameter fluctuation effects, con- 
sidering that the phase fluctuations are more important 
in underdoped cuprates due to reduced superfluid den- 
sity. Since the microscopic origin of this deformation 
is not understood, further experimental investigation on 
the temperature variation of the gap anisotrpy would be 
desirable. 

Some of the above features are shown to be shared by 
other theoretical results in the normal state counterpart. 
For instance, the form of the density of states obtained 
above is similar to that above T c when the Gaussian pair- 
ing fluctuations are incorporated.EZI Also a similar but 
much more pronounced deformation of the cZ-wave spec- 
tral gap was obtained abpiie T c using a self-consistent 
conserving approximation.E3 

The author thanks Alan Dorsey and Rob Wickham for 
helpful discussions and comments. This work was sup- 
ported by the National High Magnetic Field Laboratory 
and by NSF grant DMR 96-28926. 



13 Ch. Renner et al, Phys. Rev. Lett. 80, 149 (1998); A. K. 
Gupta and K.-W. Ng, Phys. Rev. B 58, R8901 (1998). 

14 D. S. Marshall et al, Phys. Rev. Lett. 76, 4841 (1996); 
H. Ding et al, Nature 382, 51 (1996); M. R. Norman, M. 
Randeria, H. Ding and J. C. C ampuzano, Phys. Re v. B 57, 



:ond-mat/9712100 



R11093 (1998); H. Ding et al, 

15 M. Franz and A.J. Millis, Phys. Rev. B 58, 14572 (1998). 

16 I. Kosztin, Q. Chen, B. Janko, and K. Levin, Phys. Rev. B 
58, R5936 (1998). 

17 S. Marcelja, Phys. Rev. B 1, 2351 (1970). 

18 J. R. Engelbrecht, A. Nazarenko, M. Randeria, and E. 
Dagotto, Phys. Rev. B 57, 13406 (1998). 



Z 

IT 
z 





M 
j \ 




,' \ 








\ /' \ " 








^^^^ ' / 

^ / ~~~ 


t" r 


/ / 
















.'' y — o.2 




0.5 




' 0.93 




0.98 



co/A (0) 

FIG. 1. Total density of states (N(u>)/No) as a function 
of energy at four temperatures (T/T c = 0.2, 0.5, 0.93, and 
0.98). 



* Electronic mail: hjk@phys.ufl.edu. 

1 For an overview of the experi mental results and their in- 



terpretation, see M. Randeria, cond-mat/9710223 



Miyakawa, et al, 3ond-mat/9809395 



2 See, for example, N. 
See also Ref. [12]. 

3 J. Maly, B. Janko, and K. Levin, Phys. Rev. B 59, 1354 
(1998); Y.-M. Vilk et al, J. Phys. Chem. Solids 59, 1 873 
(1998); M. Capezzali and H. Beck, |cond-mat/980621l| , to 
appear in Physica B. 

4 See also the review article, A. A. Varlamov, G. Balestrino, 
E. Milani, and D. V. Livanov, to appear in Adv. Phys. 

5 V. J. Emery an d S. A. Kivelson, N ature 374, 434 (1995); 
J. Corson et al, cond- mat/98102801 

6 H.-J. Kwon and A. T. Dorsey, |:ond-mat/9809"225| (Phys. 
Rev. B March 1999). 

7 J. R. Engelbrecht, M. Randeria, and C. A. R. Sa de Melo, 
Phys. Rev. B 55, 15153 (1997). 

8 A. A. Varlamov and V. V. Dorin, Sov. JETP 64, 1159 
(1986). 

9 R. A. Smith, M. Yu. Reizer, and J. W. Wilkins, Phys. Rev. 
B 51, 6470 (1995). 

10 M. R. Norman et al, Nature 392, 157 (1998). 

11 J. Mesot, et al, |cond- mat/9812377 . 

12 A. G. Aslamazov and A. I. Larkin, Soviet Solid State, 10, 
875 (1968). 




FIG. 2. The comparison of the density of states peak to 
the mean-field gap magnitude. The solid line is the total 
density of states peak position and the long dashed line is the 
mean-field gap. The dot-dashed line is the difference between 
the two quantities. 
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FIG. 3. The peak points of the angle-resolved density of 
states (solid line) vs the mean-field gap (dashed line) at 
T/T c = 0.93. The dot-dashed line is the DOS peak curve 
evaluated with only the order parameter amplitude fluctua- 
tions. 
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FIG. 4. The angular variation of the quasiparticle scatter- 
ing rate (r(w)/A (0) ) evaluated at w/A (T) = 1, 0.5, and 
0.3 at T/T c = 0.93. 
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